Quantum turbulence in propagating superfluid vortex front 
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We present experimental, numerical and theoretical studies of a vortex front propagating into a 
region of vortex- free flow of rotating superfluid 3 He-B. We show that the nature of the front changes 
from laminar through quasi-classical turbulent to quantum turbulent with decreasing temperature. 
Our experiment provides the first direct measurement of the dissipation rate in turbulent vortex 
dynamics of 3 He-B and demonstrates that the dissipation is temperature- and mutual friction- 
independent in the T — > limit, and is strongly suppressed when the Kelvin- wave cascade on vortex 
lines is predicted to be involved in the turbulent energy transfer to smaller length scales. 

PACS numbers: 67.57.Fg, 47.32.-y, 67.40.Vs 



Turbulent motion of fluids with low viscosity, like wa- 
ter or air is a general phenomenon in Nature and plays an 
important role in human everyday life. Nevertheless, var- 
ious features of turbulence are not yet well understood; 
thus a possibility to study turbulence from another, non- 
classical viewpoint looks promising. Flows in superfluids 
(with zero viscosity) can be turbulent; understanding this 
form of turbulence is crucial for our ability to describe 
superfluids, including practical applications like cooling 
superconductor devices. Turbulence in superfluid 4 He at 
relatively high temperatures T ~ 1 K has been studied 
for decades Fermi- fluid 3 He is very different from the 
Bose-fluid 4 He; this allows us to study in 3 He some as- 
pects of superfluid turbulence not available in 4 He Q, Q. 

In this Letter we report the first experimental obser- 
vation and study of propagating laminar and turbulent 
vortex fronts in rotating 3 He which opens a possibility to 
directly measure the rate of kinetic energy dissipation in 
a wide temperature range down to 0.18 T c of the critical 
temperature T c w 2.43 mK (at 29 bar pressure). Theo- 
retical and numerical analysis of the propagating front 
allows us to identify the mechanisms of energy dissipa- 
tion and to discover a dissipation anomaly: the satura- 
tion of the rate of turbulent kinetic energy dissipation at 
a nonzero value at temperatures below 0.22 T c . 

The idea of the experiment is as follows: in rotating 
^e (in the B-phase) we can create a vortex-free Landau 
state, which in the absence of external perturbations per- 
sists forever. In this state the superfluid component is at 
rest in the laboratory frame, while the normal component 
is in rigid rotation with angular velocity Q. The Landau 
state is met ast able, having larger free energy than the 
stable equilibrium vortex state. The latter consists of 
rigidly rotating normal and superfluid components with 
a regular array of rectilinear quantized vortices. When 
we inject a seed vortex into the Landau state, we ob- 
serve a rapid local evolution of the vorticity toward the 
equilibrium state. A boundary between the vortex- free 
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FIG. 1: Experimental setup. (Left) Configuration for vor- 
tex injection using AB interface instability in which two vor- 
tex fronts propagate independently in the upper and lower 
B phase sections (as shown by hollow arrows). (Right) Con- 
figuration for vortex injection from the orifice when a sin- 
gle vortex front travels first through the bottom and then 
through the top NMR pick-up coils. Superimposed inside the 
sample container are snapshots of vortex configurations from 
numerical simlulations of vortex expansion at two different 
temperatures. The configuration at 0.3T C displays small-scale 
structure (shown within blow-up on the right) which is absent 
at 0.5T C . 



and the vortex states propagates with constant velocity 
Vi toward the metastable region(s), see Fig. [TJ In some 
sense this phenomenon is similar to the propagation of a 
flame front in premixed fuel. The flame front propaga- 
tion can also proceed in laminar or in turbulent regimes. 
In the latter case the effective area of the front increases 
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and its propagation speed becomes higher than in the 
laminar regime. This property finds its practical use in 
combustion engines but also has been used to describe 
intensity curves of type la supernovae [4(. In all these 
cases a metastable state of matter is converted to the 
stable state in the front and Vf is determined by the rate 
of dissipation of the released energy. In our case (with 
density p s and velocity U s — Sir of the superfluid com- 
ponent in a sample container of radius R) the dissipation 
rate of total kinetic energy, £(£), is related to Vf as 

dS/dt = -7rp s VfQ 2 R 4 /4 . (1) 

Thus by measuring Vf one determines directly d£/dt in 
turbulent vortex dynamics. So far measurements in the 
turbulent regime only concerned the vortex density [2]. 

Our main experimental observation is that d£(t) / dt 
does not go to zero in the T — > limit. At high enough 
T the front is laminar and its velocity is determined 
by mutual friction between the normal and superfluid 
components Vf(T) w a(T)QR, where a(T) is a dimen- 
sionless mutual-friction coefficient. Although a(T) — > 
when T — > 0, the measured velocity Vf(T) saturates at 
a constant value which corresponds to an effective fric- 
tion a e ff ~ 0.1. We interpret this behavior to be similar 
to the viscous anomaly in classical turbulence, where the 
dissipation rate does not vanish when viscosity v goes to 
zero. The viscous anomaly appears due to cascading of 
energy to smaller length scales. When v — > the small- 
est length scale decreases but the global dissipation rate 
does not change. We believe that a similar mechanism, 
which we can call mutual friction anomaly, applies to su- 
perfluids where mutual friction plays the role of viscosity. 
One can say that the viscous and mutual friction anoma- 
lies are particular cases of a more general phenomenon, 
which we call the dissipation anomaly: The nonzero rate 
of energy dissipation in the limit of a vanishingly small 
parameter, that governs dissipation. 

As mutual friction decreases and turbulent motion 
reaches progressively smaller length scales, eventually 
quantized vortex lines become important. The energy 
cascade on length scales smaller than the intervortex 
distance and the nature of dissipation on such scales 
are currently the central questions in research on tur- 
bulence in superfluids [lj]. At the moment only theo- 
retical speculations exist on the role of the non-linear 
interaction of Kelvin waves, resulting in a Kelvin- wave 
cascade terminated by quasiparticle emission, and 

on the role of vortex reconnections which can immedi- 
ately redistribute energy over a range of scales and also 
lead to dissipation [7]. Our experiments show evidence 
for the importance of the Kelvin-wave cascade: We have 
observed a rapid decrease of the front velocity with de- 
creasing temperature in the region where the introduc- 
tion of sub-intervortex scales to energy transfer is ex- 
pected. We propose an explanation of this effect using a 
model with a bottleneck crossover between quasi-classical 



(at super-intervortex scales) and quantum turbulence (at 
sub-intervortex scales) 0] . 

Experiment: Our measurements are performed in a ro- 
tating nuclear demagnetization cryostat at Si ~ lrad/s. 
The 3 He-B sample at 29 bar pressure is contained in a 
cylindrical cell with radius R = 3 mm and length 110 mm, 
oriented parallel to the rotation axis, Fig.[TJ Pick-up coils 
of two independent NMR spectrometers near the top and 
bottom of the cell are used to monitor the vortex con- 
figuration To prepare the initial vortex- free state we 
heat the sample to about 0.75 T c (for rapid annihilation 
of all vortices) and then cool it in the vortex-free state in 
rotation to the target temperature. We can inject seed 
vortices in the middle of the cell, using the instability of 
the AB interface in rotation controlled with an ap- 
plied magnetic field. In this case two vortex fronts prop- 
agate independently up and down, arriving to the top 
and bottom pick-up coils practically simultaneously. A 
second injection technique uses remanent vortices which 
are trapped in the vicinity of the orifice on the bottom 
of the sample, Fig. 1 right. In this case the vortex front 
propagates upwards along the entire sample through both 
pick-up coils in succession. We determine the front veloc- 
ity dividing the flight distance by the flight time, as if the 
front propagates in steady-state configuration. Although 
this is not the case due to initial equilibration processes 
which follow injection, we believe that this simplification 
is justified here, since the two injection techniques for dif- 
ferent propagation lengths give the same result, as seen 
in Fig. El 

Experimental results on the front velocity are presented 
in Fig. [21 The temperature range is clearly divided in 
two regions: At T > 0.4T C the dimensionless front ve- 
locity Vf = Vf/QR ~ a. This agrees with previous 
measurements in this temperature range and can be un- 
derstood from the dynamics of a single vortex, when 
inter- vortex interactions are ignored [9]. We call this 
region the laminar regime. The measured values of Vf 
are slightly below a. We believe that the difference is 
caused by the twisted vortex state [10] behind the front. 
It reduces the energy difference across the front and cor- 
respondingly the front velocity. Estimating the reduc- 
tion factor from the uniform twist model [l(j we get 
V f ,iam = [2/ log(l + 1/q 2 ) - 2/q 2 ]a, where q = a/ (I - a') 
and a' is the reactive mutual friction coefficient. This 
dependence, (the dash-dotted line in Fig. [2]) is in good 
agreement with the experiment. 

The new feature is the behavior at T < 0.4 T c . Here Vf 
rapidly deviates to larger values than aSlR and eventu- 
ally becomes constant in the T — > limit with a peculiar 
transition from one plateau to another at around 0.25 T c . 
We attribute this behavior to turbulent dynamics and 
analyze it below in more detail. 

Simulations: To clarify the vortex front formation and 
propagation, we simulate vortex dynamics using the vor- 
tex filament model with full Bio-Savart equations and an 
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FIG. 2: Color online. Scaled front velocity vs temperature. 
Open circles and squares mark measurements with injection 
using AB interface instability and orifice trapping, respec- 
tively, at Q between 0.8 and 1.5rad/s. Dashed line is the 
measured value of mutual friction coefficient a(T) fl8l]. ex- 
trapolated below 0.35 T c with exp(-A/T) law 0. Filled 
diamonds are results of numerical simulations. Dash-dotted, 
thin solid and thick solid lines show model approximations 
which sequentially account for dissipation in the large-scale 
motion, turbulent energy transfer and bottleneck effect. (In- 
sert) Value of parameter b(T) in Eq. (|4} which was used in 
producing the thick solid line in the main panel. 



additional solution of the Laplace equation for solid wall 
boundary conditions [ll| in an ideal cylinder with length 
40 mm and diameter 3 mm rotating at Q = lrad/s. In 
the initial configuration the equilibrium number of vor- 
tices is placed close to one end plate of the sample as 
quater-loops between the end plate and the side wall of 
the sample. During evolution vortices form a front prop- 
agating along the sample (for movies see Ref. [Hj]). The 
thickness of the front is A(r) ~ r d, where d ~ 1. For 
T > 0.4T C the front is laminar with smooth vortices, 
that only twist at large scales (Fig. [TJ left). For smaller 
T the front becomes turbulent demonstrating a lot of 
small-scale vortex structure: Kelvin waves, kinks, semi- 
loops (Fig.[TJ right). These small scale fluctuations exist 
on the background of strongly polarized vortex orienta- 
tion, which is still preserved at low temperatures. The 
calculated front velocity is in good agreement with the 
measurements, see Fig. [2j In particular, the velocity is 
lower than a QR at higher temperatures and propagation 
becomes relatively much faster at lower temperatures. 

Our analytical model of the turbulent front is based 
on the quasi-classical " coarse-grained" equation averaged 
over the vortex lines 0, [l3[ , 

U + (l-a')(U- V)C/ + V/i = -rC7, T = auj eU (2) 



which describes the evolution of the superfluid velocity 
U(r,t) at scales exceeding the crossover scale i between 
hydrodynamic and kinetic regimes of turbulent motions, 
which is of the order of the intervortex separation. In 
Eq. (|2J) /i is the chemical potential and the dissipative 
term T is taken in the simplified form [lj], in which a; e f 



can be understood as an effective vorticity. 

To estimate Vf in the quasi-classical regime, described 
by Eq. ([2]), we consider the total energy dissipation in the 
front with well developed turbulence which has two con- 
tributions. The first one originates from the mutual fric- 
tion which acts on the global scale. It can be estimated 
from Eq. ([2]) as au e fK(z,r), where K(z,r) = - (u 2 ^ is 
the turbulent kinetic energy per unit mass and u is the 
turbulent velocity fluctuations (with zero mean). The 
second contribution is determined by the usual energy 
flux in classical turbulence e ~ b K 3 / 2 (r)/ L(r) at an 
outer scale of turbulence L(r). Clearly, L(r) ~ A(r), the 
thickness of the turbulent front at given radius r near the 
centerline of the cell, or, near the surface of the cell, as 
distance to surface, R — r. In the whole cell, one can use 
an interpolation formula L~ 1 (z) = A(r) _1 + (R — r) _1 . 
The natural assumption is that this energy dissipates on 
the way to small scales either due to the mutual friction 
at moderate temperatures or converges into Kelvin waves 
at the crossover scale l(r). For the classical Kolmogorov- 
41 regime b c \ ~ 0.27 [l5[. Using Eq. (pQ) we can present 
the overall energy budget as: 
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Here we accounted for the mutual-friction correction to 
the nonlinear term in Eq. ([2j) b =^> b = b(l — a') and 
used the axial symmetry to perform the integration over 
the azimuthal angle. The region with R — r < £, where 
Eq. ([2|) is not applicable, is excluded from the integration. 

In the turbulent boundary layer the kinetic energy is 
independent of the axial distance to the wall. There- 
fore qualitatively we can replace K(z,r) and aj e f(z, f) 
by their mean values across the front, K(r) and aJ e f(r). 
Also dimensional reasoning dictates A(r)cJ e f(0 — aQr 
and K(r) = c(£lr) 2 /2 with a, c ~ 1. Now Eq. ([3]) gives: 

V{ = Vf/QR ~ (2c) 3/2 6(l -a / )A + 4ca/5a, (4) 

where A = 0.2 + d[]n(R/e) - 137/60 + M/R+ ...]»1.8 
for Q = lrad/s which gives R/i ~ 17. We take b = b c \ 
and choose the parameters a = 0.2, c = 0.25 and d = 2 
to fit the measurement in the region (0.3-0.4)T c . With 
these parameters Eq. (j4]) gives Vf ~ 0.16 in the limit 
T — > (when a = a' = 0) and a very weak temperature 
dependence up to T ~ 0.45T C . However, both in the 
region of lower and higher temperatures the experiment 
shows deviation from this "plateau" (Fig. [2j). 

The reason for this deviation at T > 0.35T C is that tur- 
bulence is not well developed near the cylinder axis where 



the shear of the mean velocity, responsible for turbulence 
excitation, decreases. Therefore in the intermediate tem- 
perature region only part of the front volume is turbu- 
lent, expanding toward the axis when the temperature 
decreases. We suggest an interpolating formula between 
laminar and turbulent regimes: Vf = (V f 2 i am +V f 2 turb) 1 ^ 2 ? 
where Vf, turb is given by Eq. (|4j). This interpolation is 
shown in Fig.[2]as a thin line for T < 0.3T C and as a thick 
line for T > 0.3T C . The agreement with our experimen- 
tal observation for T > 0.3T C is good, but there is a clear 
deviation below 0.25 T c , where a < 10 -2 . The reason 
is that we did not account adequately for the quantum 
character of turbulence, which becomes important in our 
conditions at T < 0.3 T c [161 ] . which is close to the mea- 
sured transition to the lower plateau in Fig. O 

The mean free path of 3 He quasiparticles at T ~ 0.3 T c 
is close to £ while at 0.2 T c it exceeds R. This change from 
the hydrodynamic to the ballistic regime in the normal 
component may influence the mutual friction force acting 
on the individual vortices. We neglect this effect because 
for T < 0.3 T c the mutual friction is already very small 
and does not directly affect Vf. 

At these temperatures the energy flux toward small 
scales propagates up to the quantum scale £ and vortex 
discreteness and quantization effects become most im- 
portant. Even though some part of the energy is lost 
in intermittent vortex reconnect ions, the dominant part 
proceeds to cascade below the scale £ by means of nonlin- 
early interacting Kelvin waves [H, @, 0] • These waves 
are generated by both slow vortex filament motions and 
fast vortex reconnect ion events. The point is that Kelvin 
waves are much less efficient in the down-scale energy 
transfer than classical hydrodynamic turbulence which 
leads to the bottleneck effect increasing the kinetic en- 
ergy at crossover scale up to A 10 / 3 times with respect to 
its K41 value [8]. In our experiments A = \n(£/ao) ~ 10 
(ao is vortex-core radius) and the inertial interval R/£ is 
about one decade. Therefore the distortion of the energy 
spectrum due to the bottleneck reaches the outer scale, 
which leads to an essential suppression of the energy flux 
at given turbulent energy or, in other words, to a de- 
crease in the effective parameter 6, which relates e and 
K. This effect is more pronounced at low temperatures 
when mutual friction is small, thus b(T) should decrease 
with temperature. We analyze this effect with the help 
of the stationary energy balance equation for the energy 
spectrum E k in fc-space 

in which T(T), being the damping ([2]), and the en- 
ergy flux over scales e(k) are taken in the Leiht differ- 
ential approximation [17]. In the calculations we use 
L/£ = 12 as ratio of the outer and crossover scales and 
characterize the bottleneck with the boundary condition 
E k /[k 3 d(E k /k 2 )/dk] = -4 • 10 5 at the crossover scale. 



The resulting function 6(T), shown in the insert in Fig.O 
decreases from its classical value b c \ « 0.27 down to ^ 0.1 
for T < 0.2T C . Now, accounting for the temperature de- 
pendence of b(T) in Eq. (|4|), we get the temperature de- 
pendence of the quantum-turbulent front shown in Fig. [2] 
by the bold-solid green line, which is in good agreement 
with our experimental data. This significant decrease of 
the dissipation rate in the quantum regime is a conse- 
quence of the relative proximity of the outer and quan- 
tum crossover scales in our measurement. 

Conclusions: We have established that conversion of 
met ast able vortex- free rotating 3 He-B to stable state oc- 
curs via propagation of a dynamic vortex structure, a 
vortex front, whose nature depends on the magnitude 
of mutual friction dissipation. At temperatures below 
0.45 T c sustained turbulence appears in the front, pro- 
foundly affecting the vortex dynamics. Owing to the 
energy transfer in the turbulent cascade, dissipation be- 
comes temperature and mutual friction independent in 
the T — > limit. In this regime we have observed the 
influence of a quantum cascade, involving individual vor- 
tices, on the global dissipation rate. 
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